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Abstract 

Helical edge states exist in the mixed spin-singlet and -triplet phase of a noncentrosymmetric 
superconductor (NCSS) when the pair amplitude (PA) in the negative helicity band, A_, is smaller 
than the PA in the positive helicity band, A+, i.e., when the PA in the triplet component is more 
than the same in the singlet component. We numerically determine energies of these edge states as 
a function of 7 = A_/A+. The presence of these edge states is reflected in the tunneling process 
from a normal metal to an NCSS across a bias energy eV: (i) Angle resolved spin conductance (SC) 
obeying the symmetry g s (4>) = —g s (—4>) shows peaks when the bias energy equals the available 
quasiparticle edge state energy provided \eV\ < A_. (ii) The total SC, G s , is zero but modulates 
with eV for finite magnetic field H. (iii) The zero bias peaks of G s and total charge conductance, 
G c , at finite H split into two at finite eV for moderate H. (iv) At zero bias, G c and G s increase 
with H and show peaks at \H\ ~ 7-ffo where -Ho is a characteristic field. 

PACS numbers: 74.45. +c, 74.50. +r, 74.20.Rp 



I. INTRODUCTION 



Recently discovered noncentrosymmetric (NCS) superconductors such as CePt 3 Si (Ref.l) 
and Li 2 Pt 3 B (Ref.2) having strong spin-orbit interaction (SOI) among various types of un- 
conventional superconductors are of current interest in their own right. Besides, super- 
conductivity at the interfaces, such as the interface between LaA10s/SrTi03 (Ref.3), may 
also be classified as the two-dimensional NCS superconductivity due to the strong potential 
gradient. The SOI in NCS superconductors leads to mixture of the spin- singlet (s-wave) 
and triplet (p-wave) pairing 4 ; the pairing amplitude in positive (negative) helicity band is 
A + (A_) with p y — ip x symmetry. The triplet pairing occurs in both s z = —1 and +1 
channels but their chiral p-wave symmetries are conjugate 5 to each other, where s z is the 
the spin component of a triplet pair along quantization direction. Thus nonmagnetic NCS 
superconductors have potential of producing spin current without magnetic field. These 
may produce Josephson spin tunneling current 5 between two NCS superconductors and spin 
tunneling current 7 due to Andreev reflection 6 across the junction between a normal metal 
and NCS superconductor. Both the up- and the down-spin holes will be reflected in the 
Andreev process; consequently the spin polarized tunneling current flows. 

There exists helical edge mode 7 ' 8 when the superconductor has more triplet component 
with p y ± ip x symmetry than singlet component. The low energy Andreev reflection is 
mostly due to these edge modes and the incident angle dependent spin polarized current 
flows 7 through the interface. In the presence of magnetic field, the incident-angle-integrated 
current is also spin polarized. There is no helical edge mode for purely s-wave symmetry. 
The existence of zero energy Majorana fermions at the vortex state and their obeying non- 
abelian statistics 9 is also a possibility in the NCS superconductors 10 a la chiral p-wave 
superconductor 11 such as Sr 2 Ru04 (Ref.12). 

The helical edge state is present 7 ' 8 in the NCS superconductors when the magnitude of 
the triplet component of the pair amplitude is larger than the singlet component, i.e., when 
the ratio between pair amplitudes in negative and positive helicity bands, 7 = A_/A + > 
, (A_ < A + ). Applying boundary condition at the edges, Tanaka et af have found that 
the bound state energy E is proportional to transverse momentum k y for small k y . In this 
article, we numerically obtain the energy of the edge states for all permissible k y , since all 
of these have role in the tunneling process. We find that the midgap quasiparticle energy 
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(E < A_) for the edge state decreases with 7. 

Although the tunneling charge and spin conductances for purely triplet symmetry (i.e., 
for 7 = 1) have been studied by Tanaka et af, exploration for the mixed triplet and singlet 
symmetries is necessary since in the system like Li 2 Pt3B, triplet and singlet components are 
comparable 2 . We employ the method of Tanaka et af and study tunneling conductances 
for different proportionate mixture of triplet and singlet components (7 ^ 1) here in this 
article and find new and interesting consequences. The angle resolved spin current, denoted 
as g s {<t>), shows peaks at those values of incident angle for which the energy of the incident 
electron is equal to the quasiparticle bound state energy, provided the bias energy \eV\ < A_ 
and it obeys the symmetry g s {4>) = —g s {— ( t ) ) an d hence total spin conductance G s is zero at 
zero magnetic field. However, at finite magnetic field G s is finite and obeys the symmetry 
G s (eV, H) = —G s (—eV, H) = — G s (eV, —H) . The total charge conductance G c shows a dip 
at the bias energy \eV\ — A_, a zero bias peak (ZBP) at zero magnetic field, splitting of 
the peak into two at finite bias and a dip at zero bias for moderate magnetic field, and then 
reappearance of the ZBP at higher magnetic field before it eventually vanishes at very high 
magnetic field. Although G s is zero at zero magnetic field, it shows ZBP at finite magnetic 
field. The splitting of peaks and the shifting of peaks at finite bias with the increase of 
magnetic field is similar as in the case of G c . The zero bias magnitude of both G c and 
G s increases with \H\ and show peaks at \H\ ~ ^H with H = & /(n 2 £\ d ) which is the 
characteristic field where $0 is the flux quantum, £ is the coherence length and is the 
penetration length of the superconductor. 

The article is organized as follows. In Sec. II, we derive an equation for the quasiparticle 
energy of the helical edge state in noncentrosymmetric superconductor using the boundary 
condition of forming bound states. This equation is numerically solved to find the energies of 
the quasiparticle bound states. The tunneling charge and spin conductances from a normal 
metal to a NCS superconductor in absence and presence of magnetic field are formulated in 
section III. The conductances are numerically determined and the results are presented in 
Sec. IV. We summarize our results in section V. 
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II. HELICAL EDGE STATE 



We begin with the Hamiltonian for an NCS superconductor in which Cooper pairs form 
between the electrons within the same spin-split band: 

n= J2 [^Ai A c kA +(A kA 4 A c t _ kA +h. c .)] , (i) 

k,A=± 

where £ k A = 6c + Aa|k| for Rashba SOI 13 , £ k = h 2 k 2 /(2m) — /i. Here fx, m, A, k, a, and A kA 
denote chemical potential, mass of an electron, spin-split band index (±), momentum of an 
electron, coupling constant of Rashba SOI given by V so = arj k ■ cr with r/ k = yk x — xk y and 
the Pauli matrices cr, and pair potential in band A respectively. We choose ky + ik x -w&ve pair 
in both the bands, i.e., A kA = A A A k with A k = — i exp[— i0 k ]. This corresponds to triplet 
component of pair potential At = (d k - <r)ia y with d k = ^(A + + A_)ry k , i.e., the amplitude 
of the triplet component A t = |(A + + A_) and the singlet component of the pair potential 
is A s = iA s (jy with amplitude A s = |(A + — A_) (Ref.5). Therefore the superconductor 
is purely triplet with k y + ik x -wave symmetry when A + = A_, purely singlet with s-wave 
symmetry when A_ = — A + , and triplet and singlet components with equal amplitude when 
A_ = 0. Therefore the Hamiltonian (1) in the matrix form 4 read as 

H-( ^ Ak V (2) 
^ -AI k -h'_ k J 

where A k = ^ k + V so and A k = A T + A s . The solution of the Hamiltonian (2) in the bulk 
is with the energy eigen values ±a/£+ + A^_ and ±-y/|i + A 2 _, in with the Cooper pairing 
between electrons within the same spin-split band. Correspondingly, there are two Fermi 
surfaces with Fermi momenta kp = ^fma/H 2 + a/ (ma/h 2 ) 2 + 2mfi/h 2 , i.e., kp < kp. 

Consider a two-dimensional semi-infinite NCS superconductor with the edge along y- 
direction such that the edge is located at x = and the superconductor is in the region 
x > 0. We then mix two quasiparticle and two quasihole states at and near the edge. The 
corresponding wave function will have the form as 

V s (x,y) = e ik y y [e- K+x { Cl ^+e ik F* x + c 2 *P+e- ik F* x } 

+e~ K - x {d 1 ^e ik ^ x + d 2 ^e~ ik F* x }} , (3) 

where Fermi momenta along rr-direction in two spin-split bands are kp x = *J kp 2 — k 2 . 



Quasiparticle and quasihole wave functions 7 in two spin-split bands (±) are given by 
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with ^± = (E-iT + )/A + , ^ = (S-ir_)/A_, and r± = ^A 2 ± - E 2 for an edge state with 
energy E, and sin(0 ± ) = k y /kp. Here ci, c 2 , c?i, and d 2 are the corresponding weights at 
which these four quasiparticle and quasihole states mix, and k± = mT±/kp x are the inverse 
of the length scales of localized edge state for two spin-split bands. 

The boundary condition ty(x — 0, y) = determines the ratio between the coefficients a, 
b, c, and d and consequently we find an identity 
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Putting expressions of u + /v + and -u_/t>_ in Eq. (6), we find 

e 2 + a+a_ - r + r_ - iE(r + + r_) 
= \p\ [E(A_ + A + )-^(A_r + + A + r_)] 



(7) 



for positive energy quasiparticles. An equivalent equation for edge state energy is also 
derived in Ref. 7. For a purely triplet superconductor, i.e., for A + = A_, E = |/3|A + . The 
solution of Eq. (6) as a function 7 = A_/A + for (3 = 0.5,0.1 is shown in Fig. 1. The zero 
energy edge state is possible only for (3 = for all A_/A + > , (A_ < A+). There is 
no edge state for A_ = 0, i.e., when the triplet amplitude and singlet amplitude will be of 
equal magnitude. This is because the superconductivity exists only in the band of positive 
helicity as the negative helicity band becomes normal in this case. If E = A_, u_ = v_ and 
consequently f3 — ±1 which suggests |0 + | = ir/2. 

When A t < A s , the pair amplitude in the negative helicity band is negative (A_ < 0). In 
that case signs of third and fourth components of ip~ and ip^ m Eq. (5) change. Therefore 



Eq. (7) in this case reduces to 

e 2 + a + a_ - r+r_ - iE(v + + r_) 

= |/5|[-E(A + + A_)+^(A_r + + A + r_)] . (8) 

This equation does not produce any solution in the range A_ < E < — A_ except when 
the magnitudes of A + and A_ are same and the corresponding solution will be E = ±A + . 
However, these solutions do not correspond to edge state since k± = 0. Therefore, there is 
no midgap edge bound state 7 for equal or larger singlet component compared to the triplet 
component. 



III. CHARGE AND SPIN TUNNELING CONDUCTANCE 

Consider a junction between a ballistic normal (at x < 0) metal and an NCS (at 
x > 0) superconductor. The junction is characterized by an insulating barrier at 
x = with a delta-function potential V(x) = U5(x). The Hamiltonian for the nor- 
mal metal is — £k 1. In this geometry, the wave function for an electron with spin 
a (numerically ± and symbolically j or j respectively) incident from the normal metal on 
the junction is given by 

HK<r + K-°W e e- lkF * x ] (9) 

within the "Andreev approximation", where T ipl = (1, 0, 0, 0), T ^ = (0, 1, 0, 0), T tpj l = 



(0,0,1,0), T ipl l = (0,0,0,1), and kp x = y k 2 F — ky with Fermi momentum kp in the nor- 
mal metal. Here o CT)CT , a„- a and b a - a are the parallel-spin Andreev, antiparallel-spin 
Andreev, parallel-spin normal, and antiparallel-spin normal reflection coefficients respec- 
tively. The normal and Andreev reflection processes and formation of Coopar pairs inside 
the superconductor are schematically shown in Fig. 2. 

The angle resolved charge and spin tunneling conductances are thus defined to be 14 ' 15 

9c((f>) = f 1 + \ ^2 l\ a °>°\ 2 + K-^l 2 _ \ h °A 2 ~ \ b v~°\ 2 ] ) cos<f), (10) 



9s{<f>) = (^^2 a [l°^| 2 - \ a cT-a\ 2 ~ %,a\ 2 + \ba-a\ 2 ] ^ COS (11) 



respectively at zero temperature. Here the angle is defined as k y = kp sin 0. The reflection 
amplitudes can be found out by matching the wave functions and the velocity flux at x = 0: 
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(13) 



For nonzero a, the phase space of the incident electron that takes place in Andreev 
reflection gets restricted. The angles of k in two bands inside the NCS superconductor is 
restricted by — | < <f>± < |. The conservation of momentum implies kpsmcp = kp sin <f> + = 
k F sin0_. The variation of (f>± with the incident angle is shown in Fig. 3(a) for a/vp = 0.1. 
It is clear that — C < < (f> c , where C is the critical angle of incidence beyond which 
incident electron becomes totally reflected. This critical angle corresponds to + = n/2 and 
0_ = 0_ jC . The angle C decreases with the increase of a as shown in Fig. 3(b). The total 
charge and spin tunneling conductances in the unit of normal tunneling charge conductance 



G nc become 



G r .= 



g c {4) d(j) ; G s 



^nc J —S r ^nc J —a 



g a (4>) d(j). 



(14) 



We then consider the application of magnetic field H perpendicular to the plane of the 
NCS superconductor. Assuming the penetration depth is much larger than the coherence 
length of the superconductor, the corresponding vector potential in the Landau gauge may 
be approximated as A(r) = (0, — H\d exp(— x/Xd), 0) with the penetration depth A^. In a 
semiclassical approximation where the quantization of the Landau level may be neglected, 
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the quasiparticle energy becomes Doppler shifted 17 : E — > E — HA + sin (j)/H with character- 
istic filed Hq = <3>o/(7r 2 £Ad), where coherence length £ = k F / (nmA + ) (as A + is larger among 
two pair amplitudes) and $ is the flux quantum. The Zeeman coupling may be neglected 
since the energy of Doppler shift energy is very high compared to the Zeeman energy for 
large A^. In contrast, the Zeeman energy is responsible to break the degeneracy between the 
helical edge modes in quantum spin Hall systems (QSHS) 21 ~ 23 and it modulates the trans- 
port properties. The modulation of the spin conductance with H due to the Doppler shift 
in NCS superconductor as superconducting analogue 7 to the QSHS as topological system. 

In presence of small magnetic field, where the formation of Landau levels are ignored, 
the wave function in the normal side remain as superposition of plane waves as in the case 
of zero magnetic field. We also ignore the spin reflection asymmetry arising from Zeeman 
coupling in the normal side. We numerically evaluate the coefficients a's and 6's, both in the 
absence and presence of magnetic field, using Eqs.(12) and (13) and plug them into Eqs.(lO) 
and (11) to determine angle resolved charge and spin conductances. The total charge and 
spin conductances are then evaluated using Eq.(14). The numerical results are presented 
below for a fixed parameter Z = 2U/vf characterizing the effective strength of the barrier. 
However, the qualitative behavior is independent of Z as we see below. 

IV. RESULTS 

Although the NCS superconductors do not break time reversal symmetry, angle resolved 
spin conductance is nonvanishing and g s (<f>) shows peaks at those values of for which 
energy of incident electron matches with the energy of the midgap edge state. The large 
g s {4>) is due to the presence of helical edge modes 7 in NCS superconductors. We have found 
that g s {4>) depends very weakly on oi/vf around the peak position. The variation of g s (4>) 
is shown in Fig. 4 for different values of the ratio 7 = A_/A + , and two different values 
of quasiparticle energy eV for a bias voltage V across the junction. The peak in g s {4>) is 
present for \eV\ < A_. The peak shifts towards smaller |0| for larger values of 7. However, 
the total spin conductance becomes zero since g s (</>) = —gsi—'P) f° r any values of a/vp, 7, 
and eV. 

The total charge and spin conductances for different values of H and 7 are shown in 
Fig. 5. Since G c and G s are weakly dependent on a, we choose a fixed value ajvp = 0.1. 
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The charge conductance is minimum at \eV\ = A_ in the absence of magnetic field since A_ 
is the lowest energy scale in the bulk superconductor. The zero bias peak in G c at H = 
is present as is observed 16-18 in rf-wave and predicted 19,20 in p-wave superconductors. When 
the bound state quasiparticle energy E = A_, = <p c . In that case 7 = \eV ± -^sin0 c | 
in presence of bias and magnetic field. The zero bias peak remains for > but G c 
decreases with the increase of H at high magnetic field. The ZBP in G c at finite manetic 
field splits into two sharp peaks at finite bias (one at negative bias and the other at positive 
bias) and a dip in zero bias, when 7 > J^|sin0 c . The peaks shift towards higher \eV\ and 
becomes weaker on lowering \H\ so that the ZBP reappears again at a low field. Although 
the total spin conductance G s is zero at any bias, it modulates with eV at finite H. It 
has the symmetry: G s (eV,H) = —G s (eV,—H) = —G s (—eV,H). The disappearance and 
reappearance of ZBP in the magnitude of G s and the splitting of ZBP at finite magnetic 
field is similar to that of G c . 

The ZBP in G c increases initially with the magnetic field and it subsequently decreases 
creating a peak at \H\ ~ 7-f^O) i- e -j when all the midgap edge states upto the energy A_ 
take part in the conduction process. Likewise ZBP in G s also behave sameway with the 
important exception that the latter changes sign on reversing the magnetic filed direction, 
although G s is zero at H — 0. This is an extraordinary effect on the spin as well as charge 
conductances as shown in Fig. 6 by the presence of midgap helical edge states. We observe 
that the value of \H\/H at which the peaks occur decreases with 7 since A_ decreases with 
a fixed A+. For 7 = 1, the ZBP in G c is almost constant at small H (\H\ < 0.5H ) but the 
ZBP in G s changes sharply at small H as is obtained by Tanaka et al 7 . However when 7 
is small, the ZBP in both G c and G s form peaks at much smaller field. In the system like 
L^PtsB (Ref.2), spin-triplet and -singlet components are in same order which means 7 is 
small and it is estimated to be ~ 0.24. Therefore in such systems the presence of helical edge 
states will be revealed in form of peaks for zero bias charge and spin magneto-tunneling- 
conductance at as small as ~ 0.35H ~ 0.07 Tesla magnetic field for typical values of £ ~ 10 
nm and ~ 100 nm. 

In our study so far, we have chosen Z = 5 as the parameter for barrier height. Fig. 7 
shows the variation of G c at zero bias as a function of magnetic field for different values of 
Z and 7. We notice that the qualitative behaviour, in particular the positions of ZBP are 
independent of Z. The values of the tunneling conductances increase with decreasing Z, as 
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expected. 



V. SUMMARY 

To summarize, the helical edge states 7 exist in a noncentrosymmetric superconductor 
provided the triplet-pair-amplitude is larger than the singlet-pair-amplitude, i.e., when < 
7 < 1. The energies of the midgap (E < A_) edge states decrease with 7. We have 
studied the consequence of these edge states on the charge and spin tunneling conductances 
from a normal metal to a noncentrosymmetric superconductor. The angle resolved spin 
conductance g s shows peak at an angle that correspond to the conduction through the edge 
state. The g s show peaks when the bias energy \eV\ < A_. It changes sign on the reversal 
of sign of the angle since the conduction is due to helical edge states and this change of sign 
leads to zero total spin conductane G s irrespective of the bias. However, the Doppler shifted 
energy of the quasiparticles for the application of H leads to nonzero G s and it modulates 
with eV for different magnetic fields. The zero bias peak is present at high H (although G s 
vanishes at very high H). This peak splits into two (one at positive bias and the other at 
negative bias) and a dip is formed at zero bias on reduction of the field. The double peaks 
occur when 7 > JJP- sin C and they become weaker on lowering the field so that a zero bias 
peak reappear at very low field. Similarly, the disappearance and reappearance of zero bias 
peak in total charge conductance G c also occur. Moreover, G c has a dip at \eV\ = A_. 
Interestingly, the magnitude of zero bias charge and spin magneto-conductance increases 
with \H\ and form peaks at \H\ ~ 7-f/o, i-e., when all the midgap helical edge states take 
part in the conduction process. 
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FIG. 1: (Color online) The variation of edge state energy of the quasiparticles with the ratio of 
pair amplitudes betwen two spin split bands for (3 = 0.5 and 0.1. E = /3A + for A_/A + = 1 and 
E converges towards zero for all values of A_/A + . However, E = only for = 0, i.e., <f>± = 0. 
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FIG. 2: (Color online) A schematic diagram of tunneling from a normal metal to an NCS supercon- 
ductor. Up (or down) spin electrons (filled circle) incident on the junction from the normal metal 
side gets partly reflected as both spin-up and spin-down electrons as well as holes (open circle) 
in the Andreev process making Cooper pairs inside the NCS superconductor at both positive and 
negative helicity bands. A bias voltage V may be applied across the junction. 
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FIG. 3: (Color online) (a) Variation of the angles 0+ and <f>- with the angle of incidence 4> for 
otjvp = 0.1. The critical angle cp c and correspondingly the critical angle for negative helicity band, 
denoted as (j) c - are shown, (b) Variation of <p c against a/vp. 
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FIG. 4: (Color online) The variation of spin conductance g s (in the unit of G nc ) with the incident 
angle <f> for different values of 7 at quasiparticle energy E = eV = 0.1A+ (a) and 0.4A+ (b) for 
H = 0, Z = 5, and a = O.Ivf- 
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FIG. 5: (Color online) The variation of charge conductance (a-d) and spin conductance (e-h) with 
bias energy eV. The parameters Z = 5, 7 = 0.1 (a,e), 0.2 (b,f), 0.5 (c,g), and 1.0 (d,h), and 
ol/vf = 0.1 are chosen. The magnetic filed H/Hq chosen for the panels (e-h) are -0.4 (solid line), 
-0.2 (dot and long-dashed line, (dotted line), 0.2 (dashed line), and 0.4 (dot and short-dashed 
line). 15 
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FIG. 6: (Color online) (a) Charge conductance and (b) spin conductance vs. H/Hq for Z 
cx/Vf = 0.1 and at zero bias. The curves from left correspond to 7 = 1.0, 0.6, 0.2, and 0.1. 




FIG. 7: (Color online) Variation of zero bias charge conductance with H at barrier heights Z 
and 10 with 7 = 0.4 (left panel) and 0.1 (right panel) when ajVp = 0.1 
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